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Abstract: Stochastic partial differential equations (SPDEs) are the basic tool for modeling systems where noise is important. 
SPDEs are used for models of turbulence, pattern-formation, and the structural development of the universe itself. It is 
reasonably well-known that certain SPDEs can be manipulated to be equivalent to (non-quantum) field theories that nevertheless 
exhibit deep and important relationships with quantum field theory. In this paper we systematically extend these ideas: We 
set up a functional integral formalism and demonstrate how to extract all the one-loop physics for an arbitrary SPDE subject 
to arbitrary Gaussian noise. It is extremely important to realize that Gaussian noise does not imply that the field variables 
undergo Gaussian fluctuations, and that these non-quantum field theories are fully interacting. The limitation to one-loop is 
not as serious as might be supposed: Experience with quantum field theories (QFTs) has taught us that one-loop physics is 
often quite adequate to give a good description of the salient issues. The limitation to one loop does, however, offer marked 
technical advantages: Because at one loop almost any field theory can be rendered finite using zeta function technology, we 
can sidestep the complications inherent in the Martin-Siggia-Rose formalism (the SPDE analog of the BRST formalism used 
in QFT) and instead focus attention on a minimalist approach that uses only the physical fields (this "direct approach" is 
the SPDE analog of canonical quantization using physical fields.) After setting up the general formalism for the characteristic 
functional (partition function), we show how to define the effective action to all loops, and then focus on the one-loop effective 
action, and its specialization to constant fields: the effective potential. The physical interpretation of the effective action and 
effective potential for SPDEs is addressed and we show that key features carry over from QFT to the case of SPDEs. An 
important result is that the amplitude of the two-point function governing the noise acts as the loop-counting parameter and 
is the analog of Planck's constant h in this SPDE context. We derive a general expression for the one-loop effective potential 
of an arbitrary SPDE subject to translation-invariant Gaussian noise, and compare this with the one-loop potential for QFT. 

PACS: 02.50.Ey; 02.50.-r; 05.40.-Fj 



I. INTRODUCTION 



Stochastic partial differential equations (SPDEs) are an essential tool in modeling systems where noise is relevant . 
SPDEs are used for models of many macroscopic systems, from turbulence to pattern-formation H,^, to the 

structural development of the Universe itself It is known that certain SPDEs can be studied with tools 

that transform them into equivalent (stochastic) field theories which exhibit deep and important relationships with 
quantum field theory (QFT). See, for example, |l]j|,|,§ and ||l|-|l5|. 

In this paper we set up the field theoretical "minimalist formalism" for SPDEs, and demonstrate how to extract the 
one-loop physics for an arbitrary SPDE subject to additive Gaussian noise. It is important to realize that Gaussian 
noise does not imply that the field degrees of freedom undergo Gaussian fluctuations: the combined interplay of 
interactions and fluctuations will appear in the third (and higher) cumulants for the fleld 4>{x,t). Also, the limitation 
to one-loop physics is not as serious as might be supposed: Experience with quantum field theories (QFTs) has taught 
us that one-loop physics is often quite adequate to give a good description of the sahent issues ||l6|-pi|. In fact, in QFT, 
the calculation of one-loop quantities can be augmented by means of "Renormalization Group improved perturbation 
theory" , which contains most of the relevant features of the physics to all orders in the expansion parameter Jl7| , p^ . 
(This was called "magical perturbation theory" by the authors of referenc e [p^ .) Furthermore, at one- loop (and 
higher), one can also introduce the effective action and effective potential p^26|, tools that allow one to determine 
the combined effects of interactions and fluctuations on the ground state of the system. Defining and calculating 
the one-loop effective action and effective potential is straightforward. Interpreting the physical signiflcance of these 
quantities is more subtle. For arbitrary SPDEs it may not even be meaningful to deflne a notion of physical energy. 
Even when the physical energy makes sense, dissipative effects may vitiate energy conservation (even when noise is 
absent). We therefore spend some effort in establishing that certain key features of the effective action for QFTs carry 
over to SPDEs. In particular, we demonstrate that it is still meaningful to deflne and calculate the effective potential 
and look for its minima. The minima of the effective potential correspond to ground states of the system, and the 
locations of these minima are equal to stochastic expectation values of the fluctuating field in the presence of noise. 
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While it is possible to provide an abstract non-perturbative definition of the effective action |2J] , in order to proceed 
with explicit calculations (such as for the one-loop effective action) one needs a perturbative procedure based on an 
expansion in some small parameter. A well-known procedure of this type, the Martin-Siggia-Rose (MSR) formalism, 
already exists in the literature 0. The MSR formalism invokes additional (unphysical) "conjugate fields", which are 
generalizations of the fictitious fields sometimes introduced to deal with the dynamics of diffusion. These fictitious 
fields permit one to extend some of the procedures of conservative physical systems to diffusion. For instance, Morse 
and Feshbach state: "the dodge is to consider, . . ., a 'mirror-image' system with negative friction, into which the 
energy goes which is drained from the dissipative system". (See |27|] page 298.) In this paper we do not make use 
of the conjugate field formalism of MSR and, instead, proceed in a direct way in which we only have physical fields, 
(plus possibly a nontrivial functional Jacobian that can be rewritten in terms of ghost fields) . This approach simplifies 
the calculation since it halves the number of fields one has to deal with. These two alternative formalisms are very 
similar to the situation in spontaneously broken gauge field theories, where one can use two equivalent approaches 
to perturbation theory, such as "unitary gauges" versus "renormalizable gauges" : in one case the particle content is 
explicit and in the other renormalizahility is explicit. 

After setting up the path integral formalism for the characteristic functional (partition function), Z[J], we define 
both the perturbative ||2^ and non-perturbative effective action |2^. We then focus on the one-loop effective action 
and its restriction to constant (homogeneous and stationary) fields: the effective potential l2C| , |2l|j . An important 
result is that the amplitude of the noise two-point correlation function acts as the loop-counting parameter and is the 
analog of Planck's constant h in this SPDE context. 

We conclude by deriving the formula for the one-loop effective potential of a general SPDE subject to translation- 
invariant Gaussian noise. This formula has a strong resemblance to that obtained for ordinary QFTs and allows us 
to extend the use of QFT tools in the analysis of the SPDE's effective potential. We furthermore demonstrate that 
much of the physical intuition regarding the effective action in QFTs also carries over into SPDEs. Finally, we offer 
a discussion of our results. A number of more technical issues are relegated to the Appendices. 



II. STOCHASTIC PARTIAL DIFFERENTIAL EQUATIONS 

A. Elementary definitions 

Consider the class of stochastic partial differential equations of the form 

Dcp{x,t)^F[4>{x,t)]+7^{x,t), (1) 

where D is any linear differential operator, involving arbitrary time and space derivatives, which does not explicitly 
involve the field 0. Typical examples are 

Diffusion equation, (2) 
Wave equation, (3) 
Langevin equation. (4) 
The function F[(j)\ is any forcing term, generally nonlinear in the field (j). Typical examples are 

in the Kardar-Parisi-Zhang (KPZ) equation, (5) 
in reaction-diffusion-decay systems (P is a polynomial), (6) 
in "purely dissipative" SPDEs. (7) 

The forcing term will typically not contain any time derivatives, but this is not an essential part of the following 
analysis except insofar as time derivatives may complicate some of the Jacobian functional determinants that will be 
encountered below. Non-derivative terms linear in the field can be interpreted either as decay rates or (if a diffusion 
term is also present) as mass terms. They can be freely moved between the differential operator D and the forcing 
term P [(/)]. If they are considered part of the forcing term then 



D = 



D 



d_ 

dt 

V 

dt^ 



m = +2 (V0)2 

F[4>] = P[4>] 
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F[(j)] = —70 describes a decay term; 

= —i'inn?(j) describes a mass term. 



(8) 
(9) 



The function ri(x^ t) is a random function of its arguments and describes the noise that we assume is present in the 
system. For the remainder of this paper, we consider field-independent additive noise. At this stage the nature and 
probability distribution of the noise are completely arbitrary and do not need to be specified. 

The noise represents our ignorance about precise details in the dynamics of the system. It could be due, for 
example, to fluctuations intrinsic to the dynamics (as in the case of Quantum Mechanics), or it could be thought 
of as representing the dynamics of short-scale degrees of freedom which have not completely decoupled from the 
macroscopic dynamics (e.g., thermal or turbulent noise), or it could be a way of implementing ignorance of the exact 
initial or boundary conditions in the system. Noise can also be a way of summarizing the necessary truncation of 
the deterministic dynamics of a many body system when we try to describe it via a finite set of variables {e.g., a 
truncated BBGKY hierarchy). 

If we think of turning off the noise, we do not require that the non-stochastic partial differential equation D(f> — F[(j)\ 
be derivable from an action principle (i.e., the non-stochastic partial differential equation need not arise from a 
Lagrangian formalism). Nevertheless, once we include noise we demonstrate that the presence of noise automatically 
leads to a generalized action principle for the noisy system. It turns out that in the presence of Gaussian noise an 
equation of motion proportional to the factor (Dcf) — F[(j)\) can always be derived by varying a well-defined "classical" 
action, and that the solutions to this equation of motion will coincide with those of the non-stochastic equation, 
provided a certain Jacobian determinant is nonsingular (i.e., invertible). This is explained in full detail in Appendix 
B. 



B. Some typical examples 

An example of considerable interest is the reaction-diffusion-decay system where the SPDE is taken to be ||2 



^-uV^^ = P[^]+V. (10) 

This equation is used, for example, to describe the density (/)(a;, t) of some chemical species as a function of space 
and time when the chemical is subject to both diffusion (via v) and reaction or decay [via P{4'), a polynomial in the 
density field]. Expanding out the first few terms, 

P(0) = Po + Pl0 + i'20' + P30' + ---, (11) 

we can identify Pq with a constant (in space and time) source or sink, —Pi with the decay rate, and P2 with the 
reaction rate for the two-body reaction, etc. The noise accounts for random effects due to coupling to external sources, 
truncation of degrees of freedom, averaging over microscopic effects, etc. For n species of chemical reactant the field 
(j) is simply promoted to a vector in configuration space 4>{x, t) 4'i{x, t), {i = 1, . . . , n\. (The diffusion constant, 
and decay rates then become matrices, the noise a vector, and the polynomial Pi{(j)-j) a vector- valued polynomial with 
tensorial coefhcients.) 

A second well-known example is the massive KPZ equation (equivalent to the massive noisy Burgers equation) 

"T^ji|,iil 

^-yV''cj) = -ym'cj)+^{V<i>f + 1^. (12) 

In the fluid dynamical interpretation of the KPZ equation the fluid velocity is taken to be w = — V0. This model 
problem leads to a form of "turbulence" which is known in the literature as Burgulence [^ , ^. 

A third example is the enormous class of SPDEs known as "purely dissipative" systems |15[. Purely dissipative 
systems have SPDEs of the form 

These equations fall into our classification of general SPDEs as particular types of Langevin equations with D — dt 
and with a driving term that is a (functional) gradient P[0] = —SH[(l)\/d(l){x). The nomenclature "purely dissipative" 
is justified by the fact that in the absence of noise these systems satisfy 
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dH[<l>] 



(if) 



(14) 



Note that the reaction-difFusion-decay system can be interpreted as an example of a purely dissipative system if we 
take D = dt and 



H 



HDD 



[*] = 



dx. 



(15) 



On the other hand, the KPZ system is not a purely dissipative system, 

SFkpzWx) 



H{v) 



um 5{x — y) + Wx4'{x) ■ Va;(5(a; — y) 

vm^5{x - y) - X\7y4>{y) ■ \7y6{x - y) 
SFKPziHy)] 



6(p{x) 



(16) 



The class of purely dissipative SPDEs is a very wide one, but there are many SPDEs that arc not of purely dissipative 
type. We do not want to restrict attention to purely dissipative systems in this article and keep the discussion as 
general as possible. 



III. STOCHASTIC AVERAGES, CHARACTERISTIC FUNCTIONAL, FEYNMAN RULES 



We will focus on the stochastic partial differential equation 

Dq^ix,t)^ F[<j)ix,t)]+7^ix,t), (17) 

and analyze it using functional integral techniques: Feynman diagrams, the effective action, and the effective potential. 
We develop the field theory via the most direct route, with no conjugate fields present. 

We postpone to subsequent papers more technically involved approaches such as the Martin-Siggia-Rose Lagrangian 
(with its extra unphysical conjugate fields used for book-keeping purposes) and the hidden BRST supersymmetry 
implicit in these stochastic differential equations |p], p^ , p^|j34| . 

In this section we develop the necessary tools to construct the basic field theory and non-equilibrium statistical 
mechanics associated with equation (^. We will assume uniqueness of the solution to (|l^), and in order to calculate 
the characteristic functional, we will introduce an ensemble average over noise realizations, and the notion of delta 
functionals. Once the characteristic functional is available, we find it useful to introduce ghosts a la Faddeev-Popov 
before deriving the Feynman rules. We will only need to make one assumption about the noise: that it be Gaussian, 
i.e., that all its cumulants are vanishing except the first, {ri{x,t)), and second, = {ri{x,t)ri{x' ,t')). (See e.g., p5| .) 



A. Step 1: Uniqueness 

Let us assume that the partial differential equation (p7|), plus initial conditions, is a well-posed problem. Thus, 
given a particular realization of the noise, rj, the differential equation is assumed to have a unique solution which we 
designate as 

(l>soinix,t\ri). (18) 

This assumption is relatively mild but does imply that the nonlinearity is sufficiently weak so as not to drive us 
past a bifurcation point. On the other hand, it is known that noise in concert with non-linearities can lead to the 
phenomenon of delayed bifurcation in non-linear parabolic SPDEs |36| . If the partial differential equation is ill-posed, 
in the sense that the solutions are not unique, additional analysis must be developed on a case-by-case basis. A 
specific example of this behavior is spontaneous symmetry breaking in QFT, which causes the naive loop expansion 
to violate the convexity properties of the effective potential. This situation must be dealt with by an improved loop 
expansion ||37|-^. 
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B. Step 2: Ensemble average 

For any function Q{4>) of the field we introduce the ensemble average (over the noise), defined by 

(Q(0)) = / (Pyy) V[v] Qi(l)soinix,t\r,)), (19) 



where V[ri] is the probability density functional of the noise, ft is normalized to 1, but is otherwise completely 
arbitrary, that is 

I (P77) VM^l. (20) 
The symbol Vr] indicates a functional integral over all instances (or realizations) of the noise. 

C. Step 3: Delta functionals 

We next use a functional delta function to write the following identity 

4'soln{x, t\r]) = / {V(p) (j) S[(j) - 4>soln{x, t\r])\ 



(Vcb) cl) S[D(j> - F[cl)] - 77] ^fjj\ (21) 
where we have performed a change of variables and introduced the Jacobian functional determinant, defined by 

and its adjoint 

The above is just the functional analogue of a standard delta function result: If /(x) = has a unique solution at 
a; = xo, then 



a;n = \ dxx b{x - a;o) = \ dx x S{f{x)) \f{x)\ = dx x 6{f{x)) ^ f ' {x)[f ' {x)]* . (24) 



The delta function forces one to pick up only one contribution from the solution of the equation f{x) = 0, and the 
derivative is there to provide the correct measure to the integral. In the functional case the derivative becomes a 
determinant. It is in fact the Jacobian determinant associated with the change of variables from to Dcj) — It 
is now easy to see that one also has the identity 



Q{4>soin{x, t|r,)) EE j {V(f>) Q(0) S {Dcjy ~ - 77) ^/jj^ . (25) 
Furthermore, the ensemble average over the noise, equation (p^, becomes 

{Qm = j {-Dv) {-Dc^) V[7i] Q{^) S (D^ - - 7?) VjJ^. (26) 
The noise integral is easy to perform, with the result that for arbitrary stochastic averages one has 



(Q(0)) = / (P0) V[Dq^ - F[</)]] Q(0) yjjj^. (27) 



We see from this equation that the effect of the noise only appears in the stochastic average through its probability 
distribution V[D(j) — F[(j)]]. It is worthwhile to point out that the main difference, at this stage of the formalism, 



between the present "minimal" approach and that of MSR hes in the way the delta functional is handled. In MSR, 
instead of integrating directly over the noise, as is done here, the delta functional is replaced by its functional Fourier 
integral representation. This is the step wherein the conjugate field enters. If this latter route is taken, the noise 
integration can be performed exactly only for Gaussian noise. In the minimal formalism by contrast, the integration 
over the noise can be done exactly for arbitrary noise. It thus lends itself immediately for handling non-Gaussian 
systems: For general noise distributions we can explicitly write down the probability distribution for the fields as 



V[ct^]=V[D4>-F[4>]]^JJ^. (28) 

We will not further explore the possibility of arbitrary noise in this paper, since Gaussian noise (which manifestly 
does not imply Gaussian fluctuations of the fields) is already sufficiently general to be of great practical interest. 

The presence of the functional determinant is essential: it must be kept to ensure proper counting of the solutions 
to the original stochastic differential equation. In QFT this functional determinant is known as the Faddeev-Popov 
determinant and is essential in maintaining unitarity [^9|j2^, i.e., conservation of probability. In some particular cases 
the functional determinant is field- independent, and it is safe to neglect it. We discuss this more fully in Appendix A 
and in the companion papers pO| , pl| , but for the sake of generality we will carry these determinants along (with little 
extra cost) for the rest of this paper. 

D. Step 4: Characteristic functional (Partition function) 

A particularly useful quantity is the generating functional, or characteristic functional (partition function) , defined 
by taking 

Q{4>) = exp d'^x dt J{x, t) 4>{x, t) ^ (29) 
in equation (|l^). We define it as follows 

Z[J] ^exp (^j dx J{x) (j){x)^ ^ (30) 

(31) 
(32) 



j iVcf)) P[(/)]exp (^j dx J . 

j (P0) V[D(f> - F[(f\] exp (^J dx J ^Jj^, 



with an obvious condensation of notation, dx = d'^x dt. When there is no risk of confusion we will suppress the dx 
completely. This key result will enable us to calculate the effective action and the effective potential in a direct way. 

E. Step 5: Gaussian noise 

We will now make some assumptions about the noise: We assume it to be Gaussian. Without loss of generality we 
can take the noise to have zero mean, since if the mean is nonzero we can always redefine the forcing term F[(j3\ to 
make the noise have zero mean. We therefore take the noise to be Gaussian of zero mean, so that the only non-zero 
cumulant is the second order one. We do not need to make any more specific assumptions about the functional 
realization of the noise: the noise might (for instance) be white, power-law, colored, pink, l//-noise, or shot noise, 
and our considerations below apply to all of these cases. As long as the noise is Gaussian its probability distribution 
can be written as 

The characteristic functional (partition function) is thus seen from equation (p2h to be 



^ -^=1= Jm ^e.p (/ J,) exp / JiD, - F[,])G- (i.0 - F^) . 



(34) 
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This characteristic functional (partition function) contains all the physics of the model since it allows for the calculation 
of averages, correlation functions, thermodynamic variables, etc. Note that the noise has been completely eliminated 
and survives only through the explicit appearance of its two-point correlation function in the above. Since the 
characteristic functional is now given as a path integral over the physical field, all the standard machineryof statistical 
field theory (and quantum field theory) can be brought to bear [Q. See, for example ]l5| , |l9| -pl| , and — [5l| ]. 

This formula for the partition function demonstrates that (modulo Jacobian determinants) all of the physics of any 
stochastic differential equation can be extracted from a functional integral based on the "classical action" 

Sclass^cal = \j j {D4> ~ F[4>])G-\D<t> ~ F[(t>]) . (35) 

This "classical action" is a generalization of the Onsager-Machlup action The Onsager-Machlup paper dealt 

with stochastic differential equations rather than partial differential equations (mechanics rather than field theory), 
and was limited to noise that was temporally white. As their formalism was developed with the notions of linear 
response theory in mind, Onsager and Machlup assumed the "forcing term" F[(j)\ to be linear, so that both the noise 
and the field fluctuations were Gaussian. In our formalism all these assumptions can be relaxed: the forcing term can 
be nonlinear and in general the field fluctuations will not be Gaussian even if the noise is Gaussian. 



F. Step 6: Faddeev—Popov ghosts 



We mentioned previously that the Jacobian determinant is often (not always) field- independent. This is a conse- 
quence of the causal structure of the theory as embodied in the fact that we are only interested in retarded Green 
functions. The situation here is in marked contrast with that in QFT where the relativistic nature of the theory forces 
the use of Feynman Green functions {+ie prescription). As we explain in Appendix A, this change radically alters 
the behavior of the functional determinant. 

In order to avoid too many special cases, and to have a formalism that can handle both constant and field-dependent 
Jacobian factors, we exponentiate the determinant via the introduction of a pair of Faddeev-Popov ghost fields p5|,p3. 



J = det 



1 



det(27r/) 



exp 



(36) 



where / is the identity operator on spacetime. The g field is a so-called (complex) scalar ghost field. It is a field of 
anti-commuting complex variables and behaves in a manner similar to an ordinary scalar field except that there is an 
extra minus sign for each ghost loop. We also need to use the conjugate ghost field 5^ to handle the determinant. 

We should point out that if the operator D — [SF/S^)) is self-adjoint then J ^ . In this case \fjj^ reduces 
to \ J\ — \3'^\. In QFTs the relevant operators occurring in the Jacobian determinants are obtained from second 
functional derivatives of the action and are automatically self-adjoint. In contrast, for SPDEs there is no guarantee 
that D — {5F/54)) be self-adjoint, and in fact for the examples previously discussed (KPZ, reaction-diffusion-decay, 
and purely dissipative) this operator is not self adjoint. Instead we rely on the much weaker property that the operator 
D - {5F/5(i)) is real in order to write ^ J = ^fjj* = \ J\. In all cases we are interested in the relevant operators 
are not only real, but positive, so that the absolute value symbol can be ignored and the characteristic functional 
equation (^ is given by 



Z[J] = 



1 



v/dct[(27r)3G^ 



exp 



{Vcj,) iVg) iVg^) exp (^-i j j (i?0 - F[0])G-i (D0 - F [</.]) 



5' 



g exp 



J0 



(37) 



This procedure trades off the functional determinants for two extra functional integrals. The advantage of this 
procedure becomes clear when one develops the perturbation theory. (This Faddeev-Popov trick for exponentiating 
the Jacobian determinant is also essential in finding the hidden BRST supersymmetry.) 

It must be noted that neither pair of ghost field variables couple to an external source. This means they can 
only appear in internal lines in Feynman diagrams, a fact that will be used later on when we discuss loops and 
loop-counting. 

In Appendix A we take a closer look at the Jacobian functional determinant, its causal structure, and its specific 
form for local driving forces. In the latter part of this Appendix we make use of the perturbation theory based on 
Feynman diagrams to evaluate this functional determinant from another perspective. 
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G. Step 7: Feynman rules 



With the partition function in the form given above, (with two independent ghosts), it is now easy to develop a 
formal Feynman diagram expansion. We wish to treat the driving term F[(j)\ as the perturbation and expand around 
the free-field theory defined by setting F = 0. With this convention the free action is, explicitly, 

Sfree = J J ii I^-^] '^v' i^'l']} dx dy + J {i D g} dx. (38) 
There are two particle propagators in this free action, one for the (p field, and two for the ghost fields. Formally 

G,t, = [D]-\ (40) 

Here is the adjoint operator of D, defined by partial integration. (For instance if D = dt — I'V^ then = 
—dt — pV^.) In the interests of generality, we reiterate the fact that we have not assumed translation invariance for 
the noise, (although the noise is now Gaussian.) The momentum- frequency representation for the propagators is 

G^^{ki,ujr,k2,uJ2) = —jf — p -; (41) 

L>T(A:i,wi) D{k2,uJ2) 

G,tg(^,u;) = — (42) 
D{k,uj) 

The Feynman vertices come from the interaction piece of the action, which in this convention is 

= 11 {-[D4>] G-i F[cl>] + if [<^] G-i F[ct>]} dxdy- J ^^lg^ ^g^ dx. (43) 



^interaction 



The nature of the vertices (obtained by functional differentiation with respect to the fields present in the theory) 
depends on the structure of the forcing term F[(l)]. In the meantime we formally assert 

vertex: 

Gn{ki.LOi:k2,u)2) 

F[<P] - F[<P] vertex: +i ^ ^^^'^^ ; (45) 

Grj[ki,uJi;k2,u)2) 

t , . 

ghost vertex: —h 9 g- (46) 

Note that to turn these schematic Feynman rules into practical computational tools we will need to assume that 
is some specific local functional of the field 0. (Typically a polynomial or polynomial with derivatives.) When we 
define the effective action we will again see that the formalism can be successfully developed even for non-translation- 
invariant noise, and this derivation of the Feynman rules matches the generality of the definition of the effective action. 
This concludes, for now, the most general aspects of the discussion. 

When it comes to actual calculations in specific models, the majority of these models have noise that is not only 
Gaussian but is also translation invariant. In the interests of simplicity, we now (finally, and only for the rest of this 
particular section) indicate the effects of assuming translation invariance for the noise. This lets us take simple Fourier 
transforms in the difference variable x — y (more precisely: x — y and — ty) to see that in momentum-frequency 
space 

G,,(k,.) = ^ = G,{k,.)^ , (47) 

D'^{k,uj) D{k,uj) D{-k,-uj) D{k,uj) 

G,t,(fc,a;) = — (48) 
D[k,u}) 

The Feynman diagram vertices are now ((^ and g are here understood to be Fourier transformed) 
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D(k,u;) (I) F\<j>] 

F[<P] vertex: ^ ' // /"^^ (49) 



F[<j>] - vertex: 



i^MZM; (50) 



ghost vertex: ~i ff^ ^'^['^^ g. (51) 

As always, there is a certain amount of freedom in writing down the Feynman rules. It is always possible to take part 
of the quadratic piece in the total action and move it from the free action to the interaction term or vice versa. We 
have already seen that a linear term (e.g., i/'w?'4>) in the forcing function can with equal facility be reassigned to 
the differential operator D via the scheme D D — vrn^^ FW\ ^ ^[4'] ~ vrri^ct). This procedure can always be used to 
completely eliminate any linear term in F[(j)\. Similar but more complicated behaviour occurs if the forcing function 
contains both constant and quadratic pieces {e.g., a + hcjP'). With the conventions given above the interaction term 
contains (at least) a "cosmological constant", ^a^ J J G~^{x,y)dxdy, a quadratic piece, ab J J G~^{x,y)(l)^{y)dxdy, 
and a 0^ interaction. The quadratic piece could be moved into the free action at the cost of making the expression 
for the free propagator a little more complicated. This freedom in writing down the Feynman rules does not imply 
any ambiguity in the physical results: Moving quadratic pieces around from the interaction term to the free action 
will modify the Feynman rules but will not affect any physical quantities. 



IV. EFFECTIVE ACTION: LOOP EXPANSION 



In order to set up the formalism for the effective action, and its loop expansion, it is useful to first separate the 
two-point function for the noise into a shape, g2{x,y), and a constant amplitude, A, via the correspondence 



Gri{x,y)'^= Ag2{x,y). 



(52) 



For the case of Gaussian white noise this is automatically satisfied by definition: Gjj{x,y) — + A 5{x — y). For more 
general Gaussian noises (which describe, for instance, the effects of small scale degrees of freedom not fully decoupled 
from the physics of 0, such as in the case of a heat bath into which (j) has been immersed) this form of the two-point 
function allows the interpretation of the noise intensity ^ as a characterization of the bath-system coupling. This will 
become more evident when we compare the "effective potential" in noisy environments, equation ( |72| ) below, with 
the same object for zero temperature quantum field theory, equation (|7^). The normalization of the shape function 
g2{x,y) is essentially arbitrary, and any convenient normalization will suffice. 

Another advantage of singling out the intensity parameter A is that it is the loop-counting parameter for this 
formulation of SPDEs. To see this, one starts by writing the characteristic functional (with external sources rescaled 
for convenience) as 



Z[J] = 



v/dct[(27r)3G„ 



exp 



(P0) {Vg) (Pgt) exp -i 



{D<j,-F[^])g^\D<i,-F[<i,]) 
A 



9' 



g exp 



III 
A 



(53) 



The generating function (Hclmholtz free energy in statistical field theory) for connected correlation functions is defined 

by 



W[J] = +A{\nZ[J]-\nZ[Q\}. 



(54) 



The effective action (Gibbs free energy in statistical field theory) is then defined non-perturhatively in terms of W[J] 
by taking its Legendre transform [p0|,p4[ 



~W[J] 



J; 



8W[J] 
SJ 



J. 



(55) 



Here 0o is some suitable background (mean) field, which is taken to be the stochastic expectation value of (j) in the 
absence of external sources, J = 0. It is often but not always zero, and we retain it for generality. 
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The previous equation defines the non-perturbative effective action. For any specific example the previous equation 
is not very useful, and we often have to restrict ourselves to a perturbative calculation of the effective action, after 
singling out an expansion parameter. One can always develop a Feynman diagram expansion provided that the 
classical action can be separated into a "quadratic piece" and an "interacting term" , as we have already done. 

In the loop expansion the sum of all connected diagrams coupled to external sources J{x) is exactly W[J] as defined 
above, and the effective action r[0;0o] corresponds to all (amputated) one-particle irreducible graphs (IPI), that is, 
Feynman diagrams that cannot be made disconnected by cutting only one propagator. 

In the following argument, we will be considering diagrams contributing to the effective action. Recall that ghosts 
can only appear as internal lines, since they are not coupled to external sources. 

To see the role of the amplitude A as loop-counting parameter, note that each field propagator is proportional 
to A while each ghost propagator is independent of A. The vertices that do not include ghosts are proportional to 
A~^, while ghost vertices are independent of A. Thus each Feynman diagram contributing to the effective action is 
proportional to A-'^*"^*, where is the number of non-ghost propagators, and is the number of non-ghost vertices. 
But each ghost vertex is attached to exactly two ghost propagators (except for tadpole ghost loops), and each ghost 
propagator is attached to exactly two ghost vertices (except for tadpole ghost loops). In the case of tadpole ghost 
loops, exactly one propagator is attached to exactly one ghost vertex. This implies that if one assigns a factor A to 
each ghost propagator and a factor of A^^ to each ghost loop then one will not change the total number of factors 
of A assigned to the Feynman diagram. Thus the Feynman diagrams are proportional to A^^^ where / is the total 
number of (internal) propagators in the Feynman diagram and V is the total number of vertices, now including ghosts. 

It is the result of a standard topological theorem that for any graph (not just any Feynman diagram) I ~V — L—l, 
where L is the number of loops | T9| , pl| , ^ . It is then easy to see that field theories based on SPDEs exhibit exactly 
the same loop-counting properties as QFTs except that the loop-counting parameter is now the amplitude of the noise 
two-point function (instead of Planck's constant h). The only subtle part of the argument has been in dealing with 
the Faddeev-Popov ghosts, and it is important to realize that this argument is completely independent of the details 
of the differential operator D and the forcing term F[(j)]. When it comes to calculating the diagrams contributing to 
the effective action, the extra explicit factor of A inserted in the definition of W[J] above guarantees that the IPI 
graphs contribute to ipa] with a weight that is exactly A^ . This demonstrates that A is a bona fide expansion 
parameter. 

At this point, it becomes natural to make a comparison with the MSR (Martin-Siggia-Rose) formalism for the cal- 
culation of the effective action in Stochastic Field Theories, where one introduces a field conjugate to 0. Historically, 
this conjugate field first arose in setting up a variational approach to the diffusion equation, [of. Morse and Feshbach 
loc. cit.). The following remarks will help one to understand the differences and the complementarity of our approach 
to the MSR approach; the bottom line is related to technical issues associated to proving all-orders renormalizabil- 
ity ||l|,|l^,|l5| . The direct approach developed in this paper is akin to the ghost-free axial gauge of QCD or the so-called 
unitary gauge in the standard model of particle physics: This is a formalism well- adapted to isolating the physical 
degrees of freedom, at least perturbatively, but is not well-adapted to proving the all-orders renormalizability of the 
theory. (Proving one- loop renormalizability for specific theories is not too difficult, and we will address this issue in 
a pair of companion papers | |30| , |3T[ | . ) 

In analogy with the situation in QFT, one has three possible responses to this state of affairs: 

1. Use the MSR formalism for all calculations. This is comparable to using BRST-invariant versions of the standard 
model of particle physics to calculate scattering cross-sections and decay rates. (That is, overkill.) 

2. One could appeal to the fact that the SPDEs considered in this paper are hardly likely to be thought of as 
fundamental theories in the particle physics sense; These SPDEs are much more like "effective field theories", 
in that the noise and fluctuations in real physical systems are manifestations of our lack of knowledge of the 
short-distance physics. Viewed as effective theories, renormalizability is no longer the main guiding light it was 
once thought to be pC| ]. 

3. At a very practical level one can choose to be guided by experience with quantum field theories. It is well 
known that one-loop physics is often sufficient for extracting most of the physical information from a system. 
Calculations beyond one-loop, while certainly important at a fundamental level, are often more than is really 
needed. One of the great technical simplifications of one-loop physics is that, via zeta function technology, 
essentially any field theory can be regularized at one- loop without excessive complications [p3^p5[ . 



For these reasons we will now restrict our attention to a one-loop calculation (apart from the discussion of Feynman 
diagrams and the loop expansion everything up to this point has been valid non-perturbatively, while those discussions 
were still valid to all orders in perturbation theory). In the next section we calculate the one-loop effective action. 
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V. EFFECTIVE ACTION: ONE LOOP 



It is well known that the effective action for a field theory can be obtained by performing a Legendre transform on 
the logarithm of the characteristic functional (partition function). Writing 



where a is the parameter characterizing the fluctuations, one gets for the one-loop effective action (first order in a) 
r[0;0o] =5[0] -5[(/.o] + ia{lndet(52M)-lndet(52[(/.o])} + 0(a'). (57) 

Here 82'^— S'^S/S(f>{x)S(f){y) is the matrix of second order functional derivatives of the action S[(f>] (often called the 
Jacobi field operator). For QFT the loop-counting parameter a is Planck's constant ?i, and 52 is a second-order partial 
differential operator that depends on the field ((> via some potential-like term. The determinants of partial differential 
operators can be defined and calculated by a variety of techniques. The notation S[(j)o] is actually shorthand for 
5[((/)[J = 0])], and for a symmetric ground state {{(j)[J — 0]) = 0) one often has S[0] = 0. These terms contribute 
a constant offset to the effective action. In QFT these terms are interpreted as a field-independent contribution to 
the vacuum energy and are traditionally ignored, although in the context of cosmology, they contribute (sometimes 
catastrophically) to the cosmological constant. In the interest of generality we will make them explicit. When we 
consider field theories based on SPDEs, the loop-counting parameter a becomes A, which we singled out as the 
amplitude for the noise, and the bare action in equation ( p3|) is replaced by equations (^J) and (|37| ) 

S[(f\->kj J {iD(j)^F[(j)])g^\D(j)^F[(j)])}d'^xdtd''ydt'-iA {\nj + \nj^) (58) 
= Sclass^cal [<P] " iA (in ^ + In J^) , (59) 

where on the second line we have denoted by Sdassicaiiip] the double integral in the previous line. This is the quantity 
that we have previously defined as the nonlinear generalization of the Onsager-Machlup action to arbitrary Gaussian 
noise [||. 

The noise at this stage is Gaussian, and does not need to be translation invariant. We have explicitly kept the 
Jacobian functional determinant. Inserting equation ( [59| ) into the formula for the one- loop effective action [equation 
(p7|)], we obtain the following general result (applicable to any SPDE) 

r[(?!);0o] = Sclassical[4'] — Sclassical[4'o] 

+^{ilndet(52M) - ilndet(52[0o]) - ^lnj[(f>] - iln^^[0] + i\nj[(f>a] + ^liiJ^[^o]} 
+0{A''). (60) 

To make this more explicit, the fluctuation operator 52(0) {aka Jacobi field operator) is 

52^ - (d- - ^-^^) g,' (d-^-^)- {D<i> - F[0])52-^^. (61) 



Here the <— indicates that these operators should be thought of as acting to the left. Also (a;, y) is to be understood 
as a "matrix" with implicit sums over the indices x,y {i.e., integrations over the variables.) Note that if contains 
derivatives of (j) then 5F/5(j) will be a differential operator. Performing an integration by parts, this can be converted 
to a statement about the adjoint operator acting to the right, i.e., we can re-write S2[4>] as 

S2[<l>]=\D^ -TJ \92HD-Tj]- iDc^ - m)92'TTTJ- (62) 



Putting all this together gives the following one-loop result for the effective action 

r[,^;0o] = ^ j j d'^xdid'^ydt' {{D4>-F[<t>])g^\D(t>~F[(t>])]- h_A + 

Ulndet I Z?t - g-^ (d-^-1-\- {D<t> - F^g^^^ 



H J 

cj,o) + 0{A^). (63) 
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Grouping together the terms proportional to A, and using the representation of the functional determinant, enables 
us to rewrite the above in the alternative form 



r[^; M = kl j d'^xdid'^ydt' {{D(t>-F[cj,])g2\D(^-F[4>])} 



iyllndct 



I - 



SF^ 



92\D- 



SF 



{D^~F[4>])g, 




O(A^). 



(64) 



This expression for the one-loop effective action is instructive. It is made up of two contributions whose origin 
and physics are quite different. On the one hand, the first term (the generalized Onsager-Machlup term) gives a 
contribution whose form is directly related to both the noise shape factor and the non-noisy part of the equation of 
motion, including non-linearities. On the other hand, the log-determinant term is proportional to the noise amplitude 
(which we have seen is the loop expansion parameter) and its specific form depends also on the structures of D and 
as well as on properties of the noise shape function. Therefore, noise plays a central role in the physics of the 
SPDE and, as will be discussed below, particularly in the nature of the ground state of the stochastic system described 
by equation (||). 



VI. EFFECTIVE POTENTIAL: ONE LOOP 



We now concentrate on field configurations that are homogeneous and static. For such field configurations the 
effective action reduces to a quantity known as the "effective potential" . In this section we will calculate the effective 
potential, deferring the discussion of its physical interpretation (in the context of SPDEs) to the next section. 

The effective potential is defined as 



0o] 



(65) 



with (j) a homogeneous and static field configuration, and 17 the volume of spacetime. The effective potential at one 
loop is given by 



VW,M^hF^[<l^]^j d'^xdtg^^^ 



A 



SF 



A 



sf'' 



— Indet D - — - i— In dot ZJt - — 



A 

'ft 



In det 



92' D - 



6F 
50 



F[0]|y d^xdtg^^ 



5^F 



0O)+O(^2). 



(66) 



In order to turn this into a more tractable expression it is useful to introduce a frequency-momentum representation. 
First notice that 



d'^k duj 



d'^x dt g2 ^(x, ~ J (2 )'^+i '^'^^ '^^ ^ '^'^ exp[— i((jj< — k ■ x)] 
= g-\k^d,LU^O). 



(67) 



(It is clear from the formula for the one-loop effective potential equation ( |66| ) that the above integral has to be finite, 
or rendered finite by appropriate renormalizations of the noise correlation function and the parameters it contains.) 
We next make use of the following identity valid for a translation invariant operator X 



lndetX= I d'^x dt j dki duji / d''fc2 (S', ^1^1, ^1) lnX(fci, t^i) (5'*(fci, ^2) (5(cji, W2) (^2, ^^2!^', ^) 
d'^k duj 



(27r) 



-\nX{k,Lu). 



(68) 



Applying this to the one loop effective potential yields 
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duj 



■In 



i:>t(fc,w) 




D{k,Lu) - ^ 



SF 



mTTT- 



In 



(27r)''+i 
-(0-.0o) + O(^^). 

(Note that 32 plays two rather different roles above.) We now adopt the simphfying convention that 



g-\k^6,iO^0) 



J d'^x dt g2^{x,t) ^ I ^ g2^ik = 0,uj = 0). 



(69) 



(70) 



This is only a convention, not an additional restriction on the noise, since it only serves to give an absolute meaning 
to the normalization of the amplitude A. 

With these conventions, the one- loop effective potential can be written as 



60] = iF^ 



d'^kduj 



d'^k duj 



(27r)^+i 
In 



In 



D{k,Lu) - 



(27r)'i+i 




(71) 



which can be recast into 



d'^k duj 



0o) + O(^'). 



In 



~g,(k,Cu)F[^]^ 



(72) 



This formula is one of the central results of this paper. It shows that noise induced fluctuations modify the zero-loop 
piece of the potential in a way which is reminiscent of the situation in both statistical and quantum field theory. For 
example, in QFT one has |po[ , pT[| : 



d'^kduo 



(2^) 



d+l 



In 



(73) 



We see in equation ([72|) that the ground state structure of the SPDE (which we will soon see is obtained by minimizing 
V[0;(/)o]) depends on both the noise correlations and the nonlinearities induced by the forcing term. We also see 
explicitly how the noise amplitude is essential in the competition between deterministic and stochastic effects. 

The major difference between the effective potential for SPDEs and QFT, lies in the fact that for SPDEs the scalar 
propagator of QFT is replaced with a propagator which has a more complex structure for the equivalent of the "mass" 
term. This difference owes to the causal structure of SPDEs. 

Notice also that for SPDEs one can naturally adapt the noise to be both the source of fluctuations and the regulator 
to keep the Feynman diagram expansion finite. This follows immediately by inspection of (|7^) which shows that the 
(momentum and frequency dependent) noise shape function g2 will affect the momentum and frequency behavior 
of the one-loop integral. The finiteness, divergence structure, and renormalizability of this integral will very much 
depend on the functional form of 92- It is thus clear that we can use the noise shape function to regulate the integral, 
if we wish. 



VII. INTERPRETATION 



The physical interpretation of the effective action and the effective potential for SPDEs is considerably more subtle 
than that for the more usual QFTs. The situation is complicated by the fact that for a completely general SPDE 
it may not be meaningful to define a physical energy. Even when the SPDE is sufficiently special so that some 
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physical notion of energy may be defined, the system may be subject to dissipation: The physical energy need not be 
conserved, even in the absence of noise. Thus the effective action and effective potential for SPDEs are not related 
to the physical energy. This means that some of the physical intuition built up from QFTs may be misleading and it 
becomes important to reassess the notion of effective action, and effective potential to see how much survives in the 
SPDE context. 

The great virtue of the effective action and effective potential in QFT is that they contain all the information 
regarding the ground state of the system and its fluctuations: From a knowledge of the effective potential one can 
ascertain under what conditions the system will display one degree of symmetry or another. It is essential that most 
of these properties carry over to the case of SPDEs, otherwise the effective action and effective potential would be 
mere mathematical constructs without physical relevance. Fortunately the key features do in fact carry over: (1) the 
stationary points of the effective action still correspond to stochastic expectation values of the fields in the absence of 
an external current; (2) the effective potential governs the probability that the spacetime average of the field takes on 
specific values; (3) even when the notion of physical energy is lacking we will see that there is a notion of quasi-energy 
for SPDEs, with the quasi-energy being a measure of the extent to which the system has been driven away from its 
non-stochastic (zero-noise) configuration; and (4) the one-loop effective action will be demonstrated to describe the 
(approximate) probability for an initial field configuration to evolve into some final (in the asymptotic sense) field 
configuration under the influence of the stochastic noise. 



A. Equations of motion in the presence of fluctuations 



If one makes use of the definition of the effective action as a Legendre transform, it is easy to see that 

- Ji^l (74) 

where J[(/)] is that external current required in order that 

(0[J]) = 0. (75) 

In particular, by taking J = 

^^^ = ^ 0=(0[^ = O]). (76) 

Stationary points of the effective action occur at those (mean) field configurations which are zero-external-current 
stochastic expectation values of the fluctuating field. (Proof of this may be found for instance on page 65 of Wein- 
berg [po[ ). It is important to realize that one never needs to invoke the notion of energy to obtain this result. The 
QFT interpretation of this result, which we now see extends to SPDEs, is that the effective action gives the equations 
of motion once fluctuations (noise) are taken into account. (This is a non-perturbative result, not limited to the 
one-loop approximation.) 



B. Probability distribution for the spacetime average field 



We have previously seen that the probability distribution for the fluctuating fleld, considered as a function over 
spacetime, to take on the value (j){x, t) is given by the functional 

P[0] =V[D4>- F[4>\] ^JT^. (77) 
Now suppose we coarse-grain, by looking at the spacetime average of the fleld (j) as deflned by 

n^r ' 

and ask: what is the probability that this spacetime average take on a specific numerical value cjP. (For definiteness 
we impose periodic boundary conditions in space U,g and time T and interpret as the volume of the spacetime 
box. This has the technical advantage that the partition function Z[J] is then needed only for sources J that are 
strictly independent of space and time.) 
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The probability we are interested in is easily calculated to be 



Prob ( / d'^x dt , t) = (j) n^T = / (Vcj)) P[(j}] 5 [ / d'^x dt (/>(f , t) - Q.T 



{V(j}) / d\ V[(j}] exp ( i\ 
= I dX Z[J{x) ~ iX] exp(— lAilsT i 



d'^xdt (t>{x,t) - (j) risT 



(79) 
(80) 
(81) 



We now take the limit as figT becomes very large, and apply the method of stationary phase. By definition we have 

Z[J{x) = iX] = exp [nsT{iX(l3{X) - V[0(A); 0o]M}] (82) 

with the subsidiary condition 

0(A) 



iXA. 



(83) 



It is easy to demonstrate that 



Prob 



d^x dt (/)(f , t)^(j) QsT cx exp -QsT 



sis xT 



A 



O 



(84) 



Thus the effective potential governs the probability distribution of the spacetime average of the fluctuating field. 
Minima of the effective potential correspond to maxima of the probability density of the spacetime average field. The 
way we have set up the argument applies equally well to QFTs and SPDEs and makes no reference to the notion of 
physical energy. (This result is non-perturbative but approximate — it is not limited to one loop. If we take either the 
infinite volume or infinite time limits, then with probability one, the spacetime average field must equal one of the 
minima of the effective potential.) 



C. Action and quasi-energy for SPDEs 



Even though the physical energy may not be defined for arbitrary SPDEs, we nevertheless can demonstrate that 
there always exists a positive-semi-definite functional of field configurations, the tree-level action, and a related "quasi- 
energy" , whose minima correspond to maxima of the probability distribution of field configurations. 

From the way the functional formalism has been set up, we can always define and calculate an effective action and 
an effective potential even if the underlying non-stochastic version of the partial differential equation does not arise 
from a Lagrangian formulation. We have already seen that the effective action has a natural interpretation in terms of 
the equations of motion once fluctuations are taken into account, and that the effective potential governs fluctuations 
in the spacetime average of the field. We now go one step further: We distinguish two concepts of "energy", the 
"true physical energy" and the "quasi-energy" , and show that even if the physical energy is undefined (or possibly not 
useful due to dissipative effects) , the quasi-energy is still a useful measure of the extent to which fluctuations modify 
the non-stochastic equations of motion. We start from our general SPDE (|) 



and its non-stochastic version 



Dcf> = F[(t>] + r], 



Dcj) = F[ 



(85) 



(86) 



Sometimes this non-stochastic partial differential equation will arise from some Lagrangian, often it will not. 

Even if the PDE (D0 = F[(j>]) does not arise from a Lagrangian, the results of this paper demonstrate that it is 
always possible to assign a tree-level action to the stochastic system: 



S. 



classical 



i J J dxdy{D4>-F[4>])g, 



{Dcj)~F[^])>Q. 



(87) 



This classical action is positive semidefinite, and has minima (which are equal to zero) at field configurations that 
satisfy the zero-noise equations of motion. This is most obvious for white noise, when the action is a perfect square, 
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but the result is general. The noise two-point correlation function, [being an (infinite dimensional) covariance matrix] , 
is by definition positive definite. Therefore its inverse is also positive definite and similarly the (infinite dimensional) 
matrix is a positive definite operator. Thus this classical action (the generalized Onsager-Machlup action [Q) is 
always greater than or equal to zero. 

The classical action thus measures the extent to which a given field configuration fails to satisfy the zero-noise 
equations of motion; the measure of the deviation being weighted by the shape of the noise correlations. (In fact, if 
the amplitude A of the noise is set to zero, the action is identically equal to zero.) 

We now define the quasi-energy by 

S[(j)] = J Equas^[<|)] dt. (88) 

We justify calling this object the quasi-energy by the fact that if we treat it as a Hamiltonian functional, and put the 
resulting object into the partition function of an equilibrium statistical field theory, we get the generating functional 
for all the correlation functions (ignoring ghost Jacobians for the moment). Explicitly, we can write 

Equas;[4>\ = h j j d'^xd'^ydt' {D4>-F[4>])g^\Dcp-F[4>]). (89) 

Note that the quasi-energy depends both on the PDE and on the shape of the noise correlation function. If the 
amplitude A of the noise is set to zero, this quasi-energy is conserved and is exactly equal to zero. This quasi-energy 
can be thought of as a nonlinear generalization of the Onsager-Machlup "energy" to arbitrary Gaussian noiseQ. 

If we now restrict ourselves to homogeneous and static fields, and consider the effective potential as defined above, 
then by the procedures used in quantum and stochastic field theories, the effective potential (multiplied by the 
volume of space) is the stochastic expectation value of the quasi-energy {Equasi[4']) in the presence of the noise 
induced fiuctuations, and subject to the constraint (0) — (f). A proof of this result is provided on pages 72-73 of 
Weinberg | |20[ |. Though that proof is phrased in a Lorentzian-signature QFT language, it readily carries over to 
Euclidean-signature equilibrium statistical field theory. Once the physical energy is replaced by the quasi-energy, the 
proof can be extended to SPDEs as well. 

We have been able to show that minima of the effective potential also minimize the quasi-energy, and therefore the 
noise-induced deviations from the zero-noise equations of motion. 

D. Transition probabilities 

What is the probability that a certain initial field configuration (j)i (x) at time U evolves into a final field configuration 
4>f{x) at time tfl We have already developed the appropriate machinery to address this question. Indeed 



Prob(0/(f),t/;(/)i(f),tj) cx j (Vrj) Vlrj] 5[(f>soin{x,U;ri) - (j)i{x)] 5[(l)soinix,tf;ri) - 0/(f)] 

cx J {Vri) {V(f>) VM 6[(t,soin{x,t;rj) - (f>{x,t)] 6[(j>{x,U) - 0,(x)] S[(j){x,tf) - (j>f{x)] 
cx J (Vr,) {V(t>) V[v] S[D^ - F[^] - r,] 5[0(x, U) - 0,(f)] 5[0(f, tf) ~ <i>f{x)] 

cx j{Vc^)Vm-F[<i>]] \[jj'U[(j,{x,U)-U^)] mx,tf)-,Pf{x)] 

cx / {VC^) P[0] m^,U) - (i^^[x)] m^.tf) - 0/(f)] 



cx j (Vc^) exp(-5[(^]M) \[j~jU[c^{x,U) - c^,{x)] S[^{x,tf) - ^f{x)] 



CX 



[V^) exp(-5[0]M) ^JJ^. (90) 

'S,ti)=(lii{x) 



^The particular label one chooses to apply to this quantity is not important as long as one bears carefully in mind that this 
"energy" need not be the physical energy. 
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This is formally identical to the formula usually encountered in equilibrium statistical field theory, and everything so 
far is non-perturbatively correct. 

Now take a saddle point approximation: Find an interpolating field 0i„t(a;, t) that minimizes S[(jji\ and interpolates 
from (j)i{x) to (j)f{x). Perform the Gaussian integral about the saddle point. Then by definition of the one- loop 
effective action 



Prob(0/(f),i/;(^i(f),ti) « exp 



A 



(91) 



This is only a one-loop result, but it demonstrates that the effective action for SPDEs inherits many of the important 
features of the effective action for QFTs. 



E. Summary 



From the above, we see that the effective action and effective potential for SPDEs exhibit many of the key features 
of the effective action and effective potential of QFTs. This is important because it guarantees that not only is it 
relatively easy to calculate the one-loop effective potential, but also it is useful to do so: As is the case for QFTs, 
minima of the effective potential for SPDEs provide information about expectation values of the fields. The effective 
action also provides information about fluctuations in spacetime averaged fields, it gives information about the noise- 
induced deviations from the non-stochastic equations of motion, and it governs the transition probabilities whereby 
initial field configurations evolve to final field configurations. Thus, both the effective potential and the effective 
action are as useful for SPDEs as it is for QFTs. 

Furthermore, as demonstrated in recent work by Alexander and Eyink [^-psf, the effective potential is also a 
useful tool in a strong noise regime far from equilibrium. The major difference between those papers and our own 
formalism is that they work within the MSR approach. They also focus on strong noise regimes, while we emphasize 
that for many purposes a one-loop calculation is both computationally efficient and quite sufficient to extract many 
key features of the physics of the system. The two approaches are complementary, and where they overlap, they are 
in complete agreement. 



VIII. DISCUSSION 



In this paper we have developed a general and powerful formalism applicable to arbitrary SPDEs. We have shown 
how to convert arbitrary correlation functions associated with arbitrary SPDEs into functional integrals. (And for 
this first step the noise does not have to be Gaussian.) For Gaussian noise (not necessarily translation invariant) we 
have carried the formalism further, setting up the basic ingredients needed for Feynman diagram expansions with the 
noise amplitude serving as the loop-counting parameter, and defining a non-perturbative effective action in analogy 
with QFT. 

We hope to have convinced the reader that the "direct approach" developed in this paper is both useful and 
complementary to the more traditional MSR formalism |l p^ , p^ . Some questions can more profitably be asked and 
answered in this "direct" formalism. For instance, the fact that the noise amplitude is the loop-counting parameter is 
easy to establish in this "direct" formalism, but appears to have no analog result in the MSR formalism. The effective 
action gives rise naturally to the concept of an effective potential, a powerful construct well known and studied within 
the QFT context, where it serves to classify and compute ground states and allows one to investigate symmetry 
properties and patterns of symmetry breaking (both spontaneous and dynamic). An analogous construct can also 
be defined and calculated for stochastic field theories based on SPDEs, and we have done so in this work. However, 
for arbitrary SPDEs, such as those contemplated here, the notion of ground state and effective potential must be 
approached with extra care and their physical interpretation clarified. We have taken pains to do so, establishing that 
the minimum of the construct we call the effective potential corresponds to solving the full equations of motion (for 
homogeneous and static field configurations) in the presence of noise. 

We feel that the most interesting result of this analysis is a general formula for the one-loop effective potential for 
any SPDEs subject to translation-invariant Gaussian noise. This is still an extremely broad class of problems, and in 
a pair of companion papers we will specialize this analysis to two particular cases. First, we discuss the noisy Burgers 
equation (KPZ equation), where the effective potential approach immediately leads us to such interesting observations 
as the existence of dynamical symmetry breaking (DSB) and the Coleman- Weinberg mechanism pl[ . Second, we 
discuss the reaction-diffusion-decay system, and explicitly calculate the renormalized effective potential for 1,2, and 
3 spatial dimensions |^ . These are issues that are extremely difficult to address using the MSR approach. 
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APPENDIX A: JACOBIAN FUNCTIONAL DETERMINANT 



The Jacobian functional determinant is often (but not always) field independent, and can often (but not always) 
be discarded. In this Appendix we explore this issue in more detail. 



1. Causality: Retarded Green functions 



This discussion is a generalization of Rivers ||2l[| , pp. 155-156. There are also relevant comments in De Dominicis 
and Peliti jl^. Appendix B part C (pp. 370-371). See also the footnote on page 214 of Frisch and the discussion 
in Zinn- Justin (pp. 372-373 [f5|). We are interested in evaluating 



J = det{ D- 



SF 



(Al) 



To proceed we make some specific assumptions about the form of D. Let us confine attentions to the class of differential 
operators 



Qn 



Do{V). 



(A2) 



(If we take Dq — then Di is the diffusion operator while D2 is the wave operator, so this class of differential 
operators is still broad enough to cover almost everything of physical interest.) Now write 



Jn = det a" - Do 



det (5t") det 



I-Gn{ Do 



det(5r) expTrln 



SF 



SF 

I-G,ADo + — 









(1 




1] 









Here Gn is the retarded Green function corresponding to 9". Explicitly 

(t - t'V^i 
[n — 1)! 

One can easily check that this is a Green function by computing, for n > 1, 



and noting that 



(n-2)! 



(A3) 
(A4) 
(A5) 

(A6) 

(A7) 

(A8) 
(A9) 
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Finally, the retarded nature of the Green function is due to the presence of the Heaviside step function. 

The traces Tr in the formula for the Jacobian J^„ are spacetime traces. We may write this as Tr = trume tr space, 
and concentrate (for now) only on the trace over time trume- For n > 1 it is easy to see that trtime{Gn) = 0, and in 
fact that for all to > 0, trumeilGn]"^) = 0. To generalize this argument to the spacetime trace we need to make the 
assumption that F[(l){x, t)] does not explicitly contain any time derivatives. If this is the case, we can write 



5{t - 1') 



This now implies that for the spacetime trace (n > 1; to > 0) 



Tr 



Gn<Do 



H J 



= 0. 



(AlO) 



(All) 



Thus the retarded nature of the Green function, causes all the trace terms to vanish and we have the exact result 
that for n > 1 and F[(j)\ not containing time derivatives 



= det \dl' - Do 
= det(5r). 



(A12) 
(A13) 



This means that the functional determinant is simply a field independent constant. It is therefore irrelevant and 
may be discarded. (In particular, for n = 2, the stochastic wave equation, one never has to evaluate the functional 
determinant.) {Aside: This argument also works provided n > 1+ [the highest order of time derivatives occurring in 
Proving this is an easy exercise.) The partition function (characteristic functional) is now, for n > 1: 



Z[J] cx / {V(f>) exp 



■ J j{D^(^-FmG~\Dr,<l,-F[4>]) 



exp 



J0 



For n = \ the situation is almost as good. First note that 



(Gi)2(t,i')= / di Gi{t,T) GiiW) 



dt e{t-t) Q{t-t') 
{t-t') e{t-t'). 



Thus ir([Gi]^) = 0, and it is easy to show that for to > 1, (Gi) 



Gn 



that tr{[Gi 



(AM) 



(A15) 

(A16) 
(A17) 

0. The only term 



that survives is tr{Gi) = 9(0). But 6(0) is ill-defined and must be specified by some particular prescription. The 
prescription which is most useful in this context is the symmetric one wherein 0(0) is non zero and equals i. This may 
be justified by a limiting procedure as described, for example, in the text by Zinn-Justin (Chapter 4, pp. 69-70.) 
This symmetric prescription is equivalent to adopting the Stratonovich calculus for stochastic equations. Choosing 
0(0) = is equivalent to the Ito calculus. The Ito calculus simplifies the Jacobian determinant (to unity) at the cost 
of destroying equivariance under field redefinitions (th e Ito calculus explicitly breaks coordinate invariance in field 



space). See, for instance, Eyink [p6|, or Zinn-Justin |15|. We will stick with the symmetric prescription (Stratonovich 
calculus) for this paper, though suitable modifications for the Ito calculus are straightforward if at times tricky (the 
loss of reparameterization invariance under field redefinitions implies that all arguments involving a change of variables 
must be carefully re-assessed). 

Now for n — 1 only one of the trace terms in the functional determinant survives and we have, (with the assumption 
that F[(j)] contains no time derivatives) 



Tr 



Gi 



Tr 



Gi 



SF[<l>{£,t)] 

smt) 



9(0) / dt tr 

space 



9(0) Tr 



SF[<l>{x,t)] 
54>{y,t) 
5F[^{x,t)] - 



(A18) 
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This implies 



= det{dt) exp|-e(0) Tr 



Do 



sF[m] 



det{dt) exp{-irr[i:io]} exp<^-irr 



H{y) 

The first two factors in the last line are field independent and so may be discarded with the result that 

-SF[cp{x)]- 



5m 



Ji cx exp | — i Tr 
The partition function (characteristic functional) is now 

Z[J]a j{Vc^) exp j j {D^c^ - F[4>])G-\Di<i> ^ F[cl>]) 

'mm] 



exp 



■Tr 



sm) 



exp 



J<j> 



(A19) 
(A20) 
(A21) 

(A22) 

(A23) 
(A24) 



This means that for stochastic differential equations that are first order in time, the functional determinant must be 
kept. There are specific choices of the nonlinear driving term F[(j)\ that lead to even further simplifications. 

2. Jacobian functional determinant for local driving forces 

Suppose that F[(t){x)] is a local functional of the field (p. This implies that there exists a local function 
such that 



sFim.t)] 



sm,t') 

Evaluating the functional determinant now gives 



.F(0(f,t),V) S{t-t') S{x~y). 



J = exp (^-iS'^iO) J dt d'^xT 



(A25) 



(A26) 



Insofar as we trust the formal result 5'^{Q) — (see, for example, we can discard the functional determinant as 
an irrelevant constant^. 

For more general driving terms F[<f>\ one must keep the functional determinant. Nevertheless it is clear that for 
large classes of stochastic partial differential equations, including many of the most important and interesting cases, 
the Jacobian can be safely ignored. 

For differential operators D that are not of the form Z?„ discussed above, or driving forces F more complicated 
than those discussed above, one has to use other means of evaluating the functional determinant. 

For the noisy Burgers equation (KPZ system) the functional determinant can be shown to be a field independent 
constant that can be discarded. A proof of this will be provided in |3l|. For the rcaction-diffusion-decay system on 
the other hand, we find it more convenient to explicitly keep the Jacobian determinant |30|. 



^ This formal result is somewhat of a contentious issue, and we have found that it is often more useful and safer to either 
prove that the Jacobian is a field-independent constant , or to carry the Jacobian along for the whole calculation |^ . 
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3. Jacobian functional determinant via Feynman diagrams 



Let us now suppose that D = £)„, as discussed in the first section of this Appendix. Then_D(fc,a;) — (— iw)" — Do(fc) 
and the ghost propagator is 

= , . . \ (A27) 

where the retarded nature of the Green function now imphes that aU poles in the uj plane occur in the lower half of 
this plane. For each ghost loop (assuming F[(l)] contains no time derivatives) we must perform an integral over both 
frequency and momenta of the type 

jm^ f dujd^kP{k,h) 

" J (2^)d+in-i{H('^-'^i)]"-£»o(fc-fci)}' 

where all the poles (in oj) lie in the lower half-plane and the oji and ki are linear combinations of the momenta 
flowing into the ghost loop. (The function P{k, ki) is some possibly complicated function of the momenta, typically a 
polynomial, derived from SF/S(f>. There is also a set of external legs (derived from SF/5<f>) attached to each vertex of 
the ghost loop, but we do not need to know the detailed structure of these vertices to derive the expression above.) 

Since all the poles are known to lie in the lower half plane the contour of integration can be pushed out to infinity 
in the upper half- plane via the replacement w — > w + iA (A > 0), without changing the value of the integral. Thus we 
can write 

VA > • J ™ = / div d'^k Pik,h) 

J (27r)'^+in™i{H(c^-u;0-A]"-Do(A-fci)} 

Now take the limit A +oo to deduce In™ = 0. 

The only place that this argument fails is when the cj integral does not converge. This happens only for n = 1 
(first-order in time) and m = 1 (tadpole diagram) in which case we need to consider 

Xr^./ '^'''''^'^ . . (A30) 

J (27r)<«+i{-iw-£)o(fc)} 

This already reproduces the key results of the previous section: The functional determinant can be ignored for n > 1 
and for n = 1 it collapses to a single term. Performing the co integral for this remaining term we see 



/ 



Ii' = J ' l^t-^'^ - DomZZf^ (A31) 

^"^^ i {itt) (A32) 

(A33) 



(277)^^+1 

Jd'^k 

= -iS''{0). (A34) 



We conclude then, that the tadpole ghost diagram exactly reproduces the cxp{—^Tr{5F/5(f))) obtained by other 
means in the first section of this Appendix. In fact, Faddeev-Popov ghost techniques are in complete agreement with 
direct calculations of the functional determinant. This ghost-based analysis also makes clear why things are different 
in QFT. If one uses the Feynman propagator instead of the retarded propagator, there are poles on both sides of the 
real line and one cannot push the path of integration out to -l-ioo. 

APPENDIX B: EQUATIONS OF MOTION IN THE PRESENCE OF NOISE 
1. Effect of adding a small decay term 
We start with the zero-loop equations of motion for the SPDE 
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Fm = J. (Bl) 

In particular, for zero external source (J = 0) any solution of the non-stochastic bare equations of motion, D(j)—F[(j>] = 
0, is also a solution of the zero loop equations of motion. (Zero loops almost corresponds to setting the noise amplitude 
to zero and reducing the SPDE to its non-stochastic analog.) But there is a risk that the zero- loop equations may have 
more solutions than the non-stochastic bare equations. This potential problem arises if the operator — {SF^^ /S(p) is 
singular (so that it has a non-trivial null space [kernel] ) . If this operator is singular then there will be many different 
fields 4>{x, t) that correspond to a given J, making the whole Legendre transform procedure invalid. 

The best way to fix this is to add a small decay term in the system and then take the limit as the decay term 
vanishes. Specifically, take 

F[(l)] F[(j)] - £(/). (B2) 
This perturbed system has zero-loop equations given by 

g^' {D<l> ~ + £0) = J. (B3) 

Even if D"^ — (SF^/dcj)) is singular, the perturbed operator will not be, and so the perturbed equations of motion 
will have a unique solution 0soi„( J, e). It is appropriate to take the Legendre transform using this unique solution 
and consider the limit e — > at the end of the calculation. If — SF^Scf) is non-singular this does not change 
anything. If — SF/S(j)^ is singular this procedure provides a prescription for defining a unique solution to the 
zero-loop equations. 

The above complication is not peculiar to SPDEs and their associated non-quantum field theories; the same sort of 
behaviour also occurs in ordinary QFTs. For example in QED, J = corresponds to arbitrary constant electromagnetic 
field. In order to assure that J = has the unique solution F — 0, the easiest thing to do is to add a small photon 
mass. 




2. A vanishing theorem 

Suppose ~5F / 5(j)^ is non-singular, then the zero- loop equations of motion for J ~ Q are equivalent to the classical 
equations D(t)— F[(f)\ = 0. The effective action at zero-loop order, evaluated on solutions of the zero-loop equations of 
motion is exactly zero. 

Not only that, in fact, the one- loop effective action evaluated on solutions of the zero-loop equations of motion is also 
exactly zero. This happens due to the explicit occurrence of D(f) — F[(t>\ in the one-loop contribution to the effective 




action, [see (63) or (p4h], so that for solutions of the zero-loop equations of motion there is an exact cancellation 
between the Jacobian and the fluctuation operator 82- 

On the other hand, if _ j jg singular, just perturb the system with a small amount of e decay. The previous 
argument goes through for e ^ 0. (Technically as long as e is not an eigenvalue of L)t — 5F/54)^ .) Taking the limit 
e ^ justifies the extension of the vanishing result to the singular case. 

Now consider solutions of the one-loop equations of motion. These one-loop equations of motion are of the form 

g^^{Dcl>^F[4>]+ec^) = 0{A), (B4) 

with the RHS being a complicated expression. Nevertheless, we do not need to know exactly what this term is to 
deduce that evaluated at solutions of these equations of motion T[(l)soin] =0-1- 0{A^). 

This vanishing of the effective action at solutions of the one-loop equations of motion provides a useful consistency 
check on specific calculations. The underlying reason for this vanishing theorem is most easily addressed in the MSR 
formalism. In fact, it can be shown that 

r[</>;^o] = hjj {{D,ff^-F,ff[<j>]) g^^ {D,ff^~F,ff[cf\)]d''xdtd''ydt', (B5) 

with Def f and Fef / are some effective differential operator and effective driving force appropriate to the fully inter- 
acting theory. 
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